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ABSTRACT 

The presence of a small amount of hydrogen is expected in most single degen- 
erate scenarios for producing a Type la supernova (SN la). While hydrogen may 
be detected in very early high resolution optical spectra, in early radio spectra, 
and in X-ray spectra, here we examine the possibility of detecting hydrogen in 
early low resolution spectra such as those that will be obtained by proposed large 
scale searches for nearby SNe la. We find that definitive detections will require 
both very early spectra (less than 5 days after explosion) and perhaps slightly 
higher amounts of hydrogen than are currently predicted to be mixed into the 
outer layers of SNe la. Thus, the non-detection of hydrogen so far does not in and 
of itself rule out any current progenitor models. Nevertheless, very early spectra 
of SNe la will provide significant clues to the amount of hydrogen present and 
hence to the nature of the SN la progenitor system. Spectral coverage in both the 
optical and IR will be required to definitively identify hydrogen in low resolution 
spectra. 

1. Introduction 



Type la supernovae (SNe la) are now considered to be among the best "standardizable 
candles" available (Phillips et al. 1999), and they are the tool of choice for observational 
cosmology (Perlmutter et al. 1999; Riess et al. 1998; Riess 2000; Goldhaber et al. 2001). 
However, the progenitor system remains in doubt (for a review see Branch et al. 1995). The 
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most widely accepted current view is that a C+O white dwarf accretes hydrogen (or per- 
haps helium) from a companion star until it reaches the Chandrasekhar mass and explodes 
due to the thermonuclear burning of C+O. While the merging of two white dwarfs (the 
double degenerate scenario) remains a viable progenitor for SNe la, recent theoretical work 
has focused on the accretion of material via Roche lobe overflow from a MS or subgiant 
companion (hydrogen cataclysmic variables) or via wind accretion from a red-giant compan- 
ion (symbiotic stars). A recent search for radio emission appears to rule out the peculiar 
SN la 1986G as having a symbiotic star progenitor (Eck et al. 1995). Maxted et al. (2000) 
have found a candidate double degenerate system containing a 0.5 M white dwarf and the 
companion is more massive that 0.97 M Q , so that the total mass of the system exceeds the 
Chandrasekhar mass. While it is clear that double degenerate progenitors do exist, it is still 
not clear whether the merger of a double degenerate system leads to collapse (due to elec- 
tron capture, Saio & Nomoto 1998; Mochkovitch et al. 1997; Tutukov & Yungelson 1996) or 
explosion. It is also unclear whether sufficiently massive double degenerate systems exist in 
sufficient numbers to produce the observed sample of SNe la (Branch et al. 1995). Therefore 
the detection of hydrogen in a SN la would greatly expand our knowledge of the progenitor 
system. The question of the identity of the progenitor system remains a significant hurdle 
for the general acceptance of the observational cosmology results; is a gap in our knowledge 
of the binary stellar evolution; and impedes our understanding of the chemical evolution 
of galaxies. There have been attempts to detect hydrogen via looking for narrow Ha lines 
(Cumming et al. 1996), in the radio (Eck et al. 1995), and in X-rays (Schlegel & Petre 1993; 
Schlegel 1995). While all of these methods should be pursued, it seems likely that we will 
obtain large numbers of very early low resolution optical spectra of nearby SNe la from 
dedicated searches that will begin taking data soon. Thus, we examine the possiblity of 
detecting hydrogen in early SN la spectra. 



2. Methodology 

The W7 model (Nomoto et al. 1984; Thielemann et al. 1986) is a good match to normal 
SNe la spectra (Lentz et al. 2001), especially the outer layers visible during the pre-maximum 
phase of the supernova. The outer layers of W7 are unburned C+O enriched with solar 
metals, expanding at velocities, v > 15000 km s _1 , with a total mass of ~ 0.07 M . 

In the single degenerate scenario, there are two sources of hydrogen: 1) circumstellar 
matter that is accreting onto the white dwarf from the companion; 2) matter ablated from the 
companion by the supernova explosion itself (Livne, Tuchman, & Wheeler 1991; Marietta, 
Burrows, & Fryxell 2000). In the first case one expects a relatively small amount of material 
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that would most likely be detected in high resolution spectra, radio, and in X-rays. In the 
second case the amount of material is much larger and is at higher velocity. It is more 
uncertain how it would be mixed into the outer layers of the the supernova. It is the latter 
case that we examine here in a parameterized way, guided by the numerical calculations of 
Livne et al. (1991) and Marietta et al. (2000). 

We have mixed solar composition material into the unburned C+O layers at the top 
of the W7 explosion model. We have replaced up to half of this C+O mixture with solar 
composition material, while maintaining the overall density structure and luminosity. The 
models are homologously expanded to account for the expansion since explosion to epochs 
of 5, 10, and 15 days after explosion, corresponding respectively to 15, 10, and 5 days before 
maximum brightness using the canonical light curve rise time of 20 days (e.g., Riess et al. 
2000; Aldering et al. 2000). We calculate the 7-ray deposition and decay of radioactive 
isotopes, especially 56 Ni, and use the luminosities for W7 models fit to SN 1994D (Lentz 
et al. 2001) at each epoch. The spectra and 7-ray deposition are calculated using the 
general purpose radiative transfer code PHOENIX, version 11.9.0 (Hauschildt & Baron 1999). 
PHOENIX includes all of the effects of special relativity in the steady-state solution of the 
transport equation and energy balance. We have included detailed NLTE treatment for the 
ionic species: H I, He I II, C I-IV, O I III, Na I III, Mg I-IV, Al I-IV, Si I-V, S I-IV, 
Ca I-IV, Fe I- VI, and Co I III. 

3. Uniform Mixing in the Unburned Layer 

For the first set of numerical experiments, we have mixed solar composition material 
throughout the unburned, outer, C+O layers of W7. For each model, we have replaced 0.1, 
1, 10, 30, or 50% of the C+O mixture mass with an equal mass of solar composition material. 
Since the composition of the C+O material was designed to have solar-like mass fractions for 
metals heavier than oxygen the mass fraction of these heavier metals is not affected by the 
mixing process. The 0.1% and 1% models had no detectable differences from the unmixed 
models at any epoch. 

3.1. Day 15 

Our synthetic spectra for the mixing at day 15, or 5 days before maximum light, are 
shown in Figure 1. Other than some minor effects on the peaks in the blue part of the 
spectrum, the primary effect is the formation of a Ha line at the red edge of the Si II 
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absorption at 6150 A. In the comparison of synthetic spectra, we can discern the effect for 
10% and larger mixing. Outside the controlled parameters of a numerical experiment this 
effect on the spectrum of the 10% mixing model would likely not be recognized as Ha. 
The 50% mixing model does show a change in the shape of the 6150 A feature from what 
is normally seen, and would stand a reasonable chance of detection, while the 30% model 
would be marginally detectable. For the detection of hydrogen, a single spectrum at 5 days 
before maximum will not be adequate. 



3.2. Day 10 

At 10 days before maximum light, or 10 days after explosion (Figure 2), the effects 
of mixing are more pronounced and might be detected. At 10 days after explosion, the 
continuum optical depth at 5000 A, T g ^, at the base of the C+O layer, is 0.15 for the 
unmixed model, and is 0.5 for the 50% mixed model. This increase in electron scattering 
"washes out" the peaks and the troughs of the spectral features. Ha produces a distinct 
double-bottom to the Si II 6100 A feature. This is complicated by the presence of C II 
absorption at 6400 A. The C II feature is seen faintly in some SNe la, at these early epochs. 
Our previous models (Lentz et al. 2001), did not include C II in NLTE and did not show the 
feature. The 10% mixed model would likely be identified as C II in an observed SN la. The 
30% and 50% mixed models would be more likely to be identified as Ha if observed. There 
is a modest absorption feature at 4600 A for H/3 in the 30% and 50% mixed models that is 
buried in the electron scattering changes in that region of the spectrum. 



3.3. Day 5 

The earliest SNe la spectrum is for SN 1990N at 14 days before maximum light (Lei- 
bundgut et al. 1991). We have calculated hydrogen mixed models for W7 at day 5 (Figure 3), 
which is nominally at 15 days before maximum light. Again, the effects of the extra electron 
scattering from the additional electrons from the hydrogen influences the spectrum. Because 
hydrogen has one electron per proton, at a given mass density, the electron density is higher 
in the solar composition mixture, than in the original C/O layer. The effect of electron scat- 
tering is smaller at 5 days after explosion than 10 days for two reasons: the day 10 spectrum 
is bluer and the effects are concentrated in the blue; electron scattering is already strong in 
the unmixed model at day 5, t^j = 0.6 at 15000 km s -1 . The Ha feature is quite strong 
and makes a definitive double-bottomed absorption, with Si II at 6200 A, for solar mixing of 
10% and higher. The small indentation at the red peak of the Si II feature at 6000 A is due 
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to C II. The Ha feature could be identified as a very strong C II feature like what might be 
expected in a C+0 white dwarf merger scenerio, but one would also expect to see other C II 
features in that case. The formation of this feature by carbon would require a significant 
mass of carbon at high velocities, v > 15000 km s -1 , in a white dwarf merger model that 
produced otherwise normal SNe la spectra and light curves. 

Figure 4 compares the SN 1990N spectrum at day -14 (Leibundgut et al. 1991) to the 
spectra at day 5 (-15) with 50% mixing of solar composition material and without hydrogen 
mixing. Clearly the hydrogen is too red to match the feature although the extra electrons do 
help to extend the blue wing of the silicon line. C II creates an unseen bump in the middle 
of the hydrogen/silicon absorption trough. The identification of silicon mixed with carbon 
(Mazzali 2001; Fisher et al. 1996) is more probable than hydrogen. 

4. Variation of Mixing Depth 

We studied the depth of mixing of the solar material mixed into a SN la since high 
velocity material would most likely come from the companion and be mixed from the outside, 
but see Marietta et al. (2000) for a detailed discussion of the expected hydrogen distribution. 
The models in the previous section were all mixed to the depth of the transition from partially 
burned to unburned C+O material at 15000 km s -1 . In the full mixing models described 
above up to half of this material is replaced with solar composition material. It is expected 
(Marietta et al. 2000) that the amount of mixing will vary with depth and we address this 
effect here with a simple parameterization. 

We have mixed the 50% composition at days 5, 10, and 15, to depths of 25000 km s _1 
(7.2 x 10~ 4 M C+O above this velocity in W7), 20000 km s^ 1 (7.2 x lO" 3 M ), and 
15000 km s" 1 (7.2 x 10~ 2 M , as stated above). Fi gure 5 shows the region of spectra 
around the Ha feature for each of these epochs. The mixing to only 25000 km s" 1 and 
20000 km s -1 does not have detectable effects on the spectra with the minor exception of the 
20000 km s _1 model at day 5 which is marginally detectable, and shows evidence of a second 
feature affecting the Si II. An observation of an event similar to this model is unlikely to be 
conclusively recognized as evidence of hydrogen mixing in SNe la. Only the models mixed 
substantially throughout the C+O layer show clearly detectable modifications. Thomas et al. 
(2002) studied the shapes of several Si II lines near maximum light (see their Figure 6), and 
none of them showed features that would be indicative of the presence of hydrogen. 
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5. Conclusions 

The mixing of hydrogen and solar composition material into the outer layers of a Type la 
supernova model has little effect on the spectra, unless the solar material replaces a sub- 
stantial fraction of the SN la C+O ejecta. In fact the non-detection of hydrogen in SNe la 
is expected in most of the present progenitor scenarios. The effects of mixing in the outer 
layers diminish as the photosphere recedes into deeper layers of the ejecta. In our tests the 
most detectable results were obtained by replacing half of the mass of the unburned C+O in 
the outer layers with solar abundance material (for a total of 2.7 x 1CT 2 M of hydrogen). 
These models produced reasonably strong signals 5 and 10 days after the explosion, but are 
less prominent at day 15, still 5 days before maximum light. One fifth of that mass of mixed 
hydrogen is able to produce identifiable signals at 5 days after explosion and approximately 
half that mass is required 10 days after the explosion. 

The effects of additional electrons from the hydrogen atoms may also be a useful diag- 
nostic when coupled with detailed analysis of high quality data. Replacing half of the C+O 
with solar material, but to different depths, we found that only at 5 days after explosion 
did any model other than mixing to the bottom of the C+O layer provide any chance of 
detection. This is consistent with a picture based on the mass of the mixed hydrogen. Ap- 
proximately 0.02 M of solar composition material swept up from the wind of or ablated 
off of the companion is needed for a detection at 10 days after explosion (10 days before 
maximum light). The deeper it is mixed the later it will be detectable. For quantities of 
mixed hydrogen near the detection threshold, the C II feature in the red emission peak of 
Si II complicates the identification of Ha absorption. Confirmation requires several epochs 
to trace the development of the spectral features, and ideally another signal, such as narrow 
circumstellar hydrogen lines. Even earlier spectra would be advantageous. 

Additional spectral signatures that would prevent confusion of high-speed carbon ejecta 
that could mimic the Ha from mixed hydrogen are needed. Figure 6 shows the region around 
the Pa line. While the Pa line shows a pretty strong signature at the early epochs the P/3 line 
is not easily discerned, thus, combined IR and optical data would be likely be required for a 
definitive detection. Bowers et al. (1997) and Hernandez et al. (2000) compared the IR data 
for SNe la and most of the data is obtained at late times (SN 1998bu has data at about 13 
days after maximum); the data are obtained much later than we desire and without detailed 
spectral analysis it would be difficult to draw any conclusions. SN 1999ee (Hamuy et al. 
2002) does have very early IR spectra (approximately 11 days before maximum). It shows 
no sign of any feature near P/3, and the region around Pa is not covered spectroscopically. 
Nevertheless, early IR spectra are obtainable and should be vigorously persued. 

The results of Marietta et al. (2000) suggest that only about 10~ 4 M Q of hydrogen would 
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be stripped from the companion at v > 15, 000 km s _1 , however, they find that 10~ 3 M of 
hydrogen would be stripped from the companion at v > 10,000 km s -1 . Thus, our results 
indicate that very early photospheric spectra of SNe la might be able to detect the presence 
of hydrogen, and certainly could be used as a signal to activate high resolution, radio, and 
X-ray programs. 

Continued work in the radio, X-ray, and high resolution optical spectra can verify and 
reinforce any putative hydrogen detections. 
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Fig. 1. — The fully mixed models at Day 15. The legend labels the percentage of hydrogen 
mixing which takes the values (0,0.1,1,10,30,50)% in this and the next two figures. 




Fig. 2. — The fully mixed models at Day 10. The legend labels the percentage of hydrogen 
mixing. 
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Fig. 3. — The fully mixed models at Day 5. The legend labels the percentage of hydrogen 
mixing. 
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Fig. 4. — The SN 1990N spectrum at day -14 (Leibundgut et al. 1991) is compared to the 
spectra at day 5 (-15) with and without hydrogen mixing. Clearly the hydrogen is too red to 
match the feature, although the extra electrons do help somewhat to extend the blue wing 
of the silicon line. The identification of high velocity carbon (Fisher et al. 1996) is more 
probable than hydrogen. 
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Fig. 5. — The effects of mixing depth variation on the region around the Ha line at days 5, 
10, and 15. Half of the C/O layer has been replaced with solar composition material which 
has been mixed to depths of (15000, 20000, 25000) km s -1 as described in the legend. 
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Fig. 6. — The effects of mixing depth variation on the region around the Pa line at days 5, 
10, and 15, for various mixing fractions (mixed to 15,000 km s~ x ) as described in the legend. 



